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Abstract: 

Recently data from PAMELA, ATIC, FERMLLAT and HESS show that there 
are excesses in the cosmic ray energy spectrum. PAMELA observed excesses 
only in e+, but not in anti-proton spectrum. ATIC, FERMI-LAT and HESS 
observed excesses in e+ + e~ spectrum, but the detailed shapes are different 
which requires future experimental observations to pin down the correct data 
set. Nevertheless a lot of efforts have been made to explain the observed 
excesses, and also why PAMELA only observed excesses in e+ but not in anti- 
proton. In this brief review we discuss one of the most popular mechanisms 
to explain the data, the dark matter annihilation. It has long been known 
that about 23% of our universe is made of relic dark matter. If the relic dark 
matter was thermally produced, the annihilation rate is constrained resulting 
in the need of a large boost factor to explain the data. We will discuss in 
detail how a large boost factor can be obtained by the Sommerfeld and Briet- 
Wigner enhancement mechanisms. Some implications for particle physics model 
buildings will also be discussed. 

Keywords: Dark matter, excess, Annihilation, Boost factor, Particle 



1. Introduction 

Recently several experiments have reported excesses in cosmic ray energy spec- 
trum. Last year the PAMELA collaboration reported 6+ excesses in the cosmic 
ray energy spectrum from 10 to 100 GeV, but observed no anti-proton excess [1,2] 
compared with predictions from cosmic ray physics [3-6]. These results are com- 
patible with the previous HEAT and AMSOl experiments (e.g., Ref. [7,8]) but with 
higher precision. Shortly after the ATIC and PPB-BETS balloon experiments have 
reported excesses in the e+ -I- e~ spectrum between 300 and 800 GeV [9, 10]. The 
ATIC data show a sharp falling in the energy spectrum around 600 GeV. Newly 
published result from FERMI-LAT collaboration also shows excesses in the e+ + 
energy spectrum above the background [11]. However, the spectrum is softer than 
that from ATIC. In addition, the HESS collaboration has inferred a flat but statisti- 
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Fig. 1. Observational data and background estimate on e"*" (left) and p (right) energy spectra. 
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Fig. 2. Observational data and background estimate on e"*" -|- e energy spectrum. J is the total 
flux of + e~ . 

cally limited e+ -I- e~ spectrum between 340 GeV and 1 TeV [12] which falls steeply 
above 1 TeV [13]. The summary of data are shown in Fig. [1] and Fig. [5] adapted 
from references [6] and [11]. 

Although astrophysics calculations of background spectrum in our galaxy 
have errors due to model parameters [3-6], within reasonable ranges it is not pos- 
sible to eliminate the excesses in the energy range from 10 GeV to 1 TeV. The 
excesses in cosmic ray have generated much excitement in particle physics commu- 
nity because dark matter (DM) , which contributes about 23% energy density of our 
universe with properties different than those of the standard model (SM) particles, 
can provide a nature explanation [14-24] [25] [26-30] [31-48] [49-51] [52-54] [55-78] 
[79-95] [96-130] [131-133] [134-137]. 

If the data from recent PAMELA, ATIC, FERMLLAT and HESS are confirmed, 
one can extract a lot of information about dark matter. The mass of the annihilating 
DM serves as the cut off scale of the spectrum, the lepton spectra must have 
a cut off energy at the DM mass itid. The FERMI-LAT and HESS data would 
require that the DM mass to be around 1 to 2 TeV. The DM belongs to the weakly 
interacting massive particle (WIMP) category. 
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To produce large enough excesses with the annihilation mechanism, it requires 
modifications of the usual DM properties. This is because that for the usual DM the 
annihilation rate producing the excess signal is also related to the annihilation 
rate producing the cosmological relic DM density. The latter requires that the ther- 
mally averaged annihilation rate < cru > to be 3 x 10~^^ cm'^s"^. This annihilation 
rate is too small by a factor of 100 to 1000 to explain the observed excess. There is 
the need to boost up the spectrum with a boost factor [4-6] B of 100 to 1000. 

Several mechanisms have been proposed to produce a large boost factor, in- 
cluding the DM substructures [134-137], non-thermal production of DM mech- 
anism [131-133], the Sommerfeld mechanism [15] [25] [138-143], and the Breit- 
Wigner mechanism [15] [30] [144,145]. There are also some other proposals [146- 
149]. Detailed calculations based on the N-body simulation show that the boost 
factor from DM substructures is generally less than [134-137] ^ 10. The reason 
why non-thermal production of DM in the early universe can explain excesses is 
because that in this scenario the interaction rate responsible for the excesses are not 
directly related to the relic DM density. One basically assumes that the interaction 
rate is that required from the excess data, and therefore there is no need of a 
boost factor. The Sommerfeld and Breit-Wigner mechanisms are more of particle 
physics answers to the problem which require the existence of new particles. In order 
for the Sommerfeld mechanism to be in effective, the new particle needs to be light 
to allow long range interaction between DM. The Breit-Wigner mechanism works 
if the annihilation of DM is through s-channel and the new particle to have a mass 
twice of the DM mass. This is a resonant effect. 

If the DM is not stable and decay predominately into leptons on time scale 
longer than the age of the universe, DM decay can also provide another alternative 
explanation to the data [6] [96-130]. The scale of the mass then provides a natural 
cut-off scale mi)/ 2 for energy spectrum. To explain FERMI-LAT and HESS data, 
the DM mass is then required to be around 3 TeV. The typical life-time required to 
fit data is a few times of lO^^s. This time is much longer than the life-time of the 
universe and will not cause other cosmological problems. 

One, of course, should not exclude the possibility that there may be other ex- 
planations [150-163] One of such possibilities is that the excess is produced by 
near by pulsars. Electrons in the intense rotating magnetic field that surrounds the 
neutron star can emit synchrotron radiation that is energetic enough to produce 
electron and positron pairs, but much harder to produce proton and anti-proton 
pairs. This provides a nice answer to why PAMELA only observed positron but not 
anti-proton excesses. The resulting positron spectrum can be modeled as a product 
of a power law and a decaying exponential with a cut-off in energy. This explains 
why the spectrum falls off at higher energies. 

In this brief review, we will concentrate on discussions of how DM annihilation 
can explain excess in cosmic ray energy spectrum. The review is arrange as the 
following. In Sec. II, we review propagation mechanism relating the sources of and 
the detected spectrum. In Sec. III. we discuss excesses from DM annihilation. 
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Both Sommerfeld and Briet-Wigner mechanisms for a large boost factor will be 
explained in some details. In Sec. IV some model building aspects of DM will be 
discussed. Finally in Sec. V we give our conclusions. 

2. The cosmic spectrum and the boost factor 

The detected spectrum of on the earth are different than that of the spectrum 
produced at the sources. The propagation from the sources to the point of detection 
will distort the shape because charged particles propagate diffusively in the galaxy. 
Some of the main effects affecting the spectrum arc the interactions with interstellar 
media when going through galactic turbulent magnetic field and radiation which 
lead to energy losses of the propagating particles, and the overall convection driven 
by the galactic wind and re-acceleration due to the interstellar shock waves. It 
is a non-trivial matter to get reliable estimate for energy spectrum taking all 
effects into account. Nevertheless theoretical efforts have been made to estimate the 
background and DM signal. A commonly used numerical package is the GALPROP 
which takes into account many known astrophysics effects of our galaxy [3]. This 
package can provide many details for energy spectra of e^, anti-proton and other 
particles. To have some understanding of the physics involved and also to have some 
simple estimate of the spectrum, we will describe a simplified method to evaluate 
the spectrum from DM annihilation in the following. 

Neglecting convection and re-acceleration effects, the flux per unit energy of 
ultra-relativistic positron or electron is given by ^e(t,r,E) = f{t,r,E)/A'K with / 
obeys the diffusion equation [14] 

where K{E) is the diffusion coefficient which are usually parameterized as 
KQ{E/GeYY, and b{E) = E"^ / {GeVrE) is the energy loss coefHcient with te = 
IQ-'^^s. These terms describe transport through the turbulent magnetic fields and 
energy loss due to synchrotron radiation and inverse Compton scattering on galac- 
tic photons. Qe is the source term given by 

where dN^/dE is the spectrum of the electron or positron produced by DM anni- 
hilation, and p{r) is DM density at r. 

The above equation is then solved in a diffusive region with the shape of a solid 
cylinder that sandwiches the galactic plane, with height 2L in the z direction and 
radius R = 20 kpc in the r direction. The location of the solar system corresponds 
to r = [rsun] Zsun) = ((8-5 ± 0.5)fcj3c; 0). The boundary conditions are usually set 
to be that the or e~ density vanishes on the surface of the diffusive cylinder, 
outside of which turbulent magnetic fields can be neglected so that positrons freely 
propagate and escape. 
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The values of the propagation parameters 5, Kq and L are deduced from a variety 
of cosmic ray data and modclizations. The following three sets of parameters have 
been used frequently [164]: i) Min : 5 = 0.55, Kq = 0.000595kpcV^/y?-, L = Ikpc ; 

> = 0.46, 



Med 



0.70, Kq = OmUkpcyMyr, L = 4kpc ; and iii) Max 



Ka = 0.0765kpc7Myr, L = 15kpc . 

To finally obtain the solution for eq.(Il|) one has also to know the details of the 
DM halo profile. There are different models. Some of the popular ones can be casted 
into the following form: 



p{r) 



^ ( 1 + {r sun I TsT 



(/3-7)/a 



(3) 



where psun is the DM density at the earth position which is believed to be in the 
range 0.2 ~ 0.7 and most of the studies use 0.3. For the other parameters a, /?, 
7 and in the profile, different values have been used. For example: a) the Core 
Isothermal (CI) model [165] has = 2, /3 = 2, 7 = and = 5A:pc; b) the Navarro, 
Frenk and White (NFW) model [166] has a = 1, /3 = 3, 7 = 1 and r,, = 2Qkpc; and 
c) the Moore model [167] has a = 1, /3 = 3, 7 = 1.16 and Vg = 30fcpc. 

For practical uses, it is convenient to factor out the galactic astrophysics mod- 
eling for the propagation and particle physics modeling of the inject spectrum 
dNe/dE, and write the final flux of stationary solution of eq.ll]) [df/dE = 0) 
at the detection point in the following form 



^e{E,rsun) = B 



1 1 



47r6(£') 2\mD 



dE' <av>%^Ia{\D{E,E')) , (4) 



E 



dE' 



where A|, 



4,KoTEiiE'/GeYY-^ - (S/GeV)''-i)/((5 - 1). The function encodes 
the galactic astrophysics. The particle physics producing the inject positron is con- 
tained in dNe/dE. Numerical solutions for la have been obtained in Ref. [14] using 
the CI, NFW and More models for the Min, Med and Max parameter sets. The 
factor B is the so called boost factor. If the model of propagation is correct, the 
factor B should be equal to 1. Approximate analytic forms for have also been 
obtained in the literature. For example, for the NFW profile, the annihilation flux, 
to a good approximation one can write la in the following form [14] 



^i^o) = flo + aitanh ( — ) [ a2exp 



Cl 



C2 



as 



(5) 



where I = logio{XD/kpc). 

Fitting numerical results, the following are obtained in Ref. [14] 



Min : a = -0.9716 , b = -10.012 ; 

ao = 0.500 , ai = 0.774 , as = -0.448 , ag = 0.649 
bi = 0.096 , 62 = 192.8 , a = 0.211 , C2 = 33.88 . 



(6) 
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Med : a = -1.0203 , b = -1.4493 ; 

ao = 0.502 , ai = 0.621 , aa = 0.688 , 03 = 0.806 , 

61 = 0.891 , 62 = 0.721 , ci = 0.143 , C2 = 0.071 . (7) 

Max : a = -0.9809, 6 = -1.1456 ; 

ao 0.502 , ai = 0.756 , 03 1.533 , 03 ^ 0.672 , 

bi ^ 1.205 , 62 0.799 , ci ^ 0.155 , C2 = 0.067 . (8) 

To compare with data, one has also to have knowledge about the background. 
The background fluxes from astrophysical sources are believed to be mainly 
due to supernova explosions for the primary electrons and from the interactions 
between the cosmic ray nuclei, such as proton and light atoms, such as hydrogen 
and helium, in the interstellar medium for the secondary electrons and positrons. 
They are commonly parameterized in the following form [5], 

0.16^-11 0.7^"-^ 



^bkgd.prini U.IUJIV ^bkgd.sec 



1 + ll£:0-9 + 3.2£;2 i5 ' " 1 + llO^i-s + 580£;4-2 ' 

(pbkgd.sec _ (q) 

1 + 650£;2.3 + i5000£;4.2 • V ) 

In the above the energy E is in unit GeV. 

With the background spectra known, one can say more about the role of DM 
in explaining the observed excesses. 

The above set of background spectra agree well with the more sophisticated 
numerical simulation results shown in Figs. [1] and [2] for energy larger than 10 
GeV. The e± excesses observed by PAMELA, ATIC, FERMLLAT and HESS have 
to be due to other contributions. If DM annihilation can explain the excesses, 
the ratio <I>g+ / ($6+ + *&e- )data observed by PAMELA, and the normalized flux 
E^{^^+ + <^>e-)data obscrvcd by ATIC, FERMLLAT and HESS must be equal to 
the background plus the DM generated fluxes. 



^ \ ^ ^bkgd^sec 



(10) 



$„+ + $„_ / , , $D _L ^bkgd.sec , $n , ^bkgd,sec ^^bkgd,przm 

+ $e-)data = Se^($f+ + ^^^^d.-c ^ ^ ^bkgd^sec ^ ^^bkgd,pr^m~^ 

Note that in the above a parameter k has been introduced in the equations to take 
care uncertainties of primary background e~ production. By adjusting k one expects 
to make a better fit to simulation data and is usually taken to be [168, 169] 0.8. If 
there is no DM contribution, the data show a large excesses. DM annihilation is 
one of the most interesting possibilities. If true the excess data then determine 
how DM contribute to the cosmic spectrum through dNe/dE. 

For DM annihilation, the parameters involved regarding DM properties are the 
DM mass m/j, the annihilation rate < av >, and the spectrum dNf./dE from 
annihilation. Since the excesses go up to the TeV region, mo must also be in the TeV 
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region to cover the full range of excesses. With DM mass fixed, one needs to worry 
about what type of final states in the DM annihilation can produce the observed 

spectrum shape. This depends on the properties of DM that is responsible for 
dNe/dE. Several model independent studies have been carried out [15-24] and find 
that the shape of the spectrum can be easily obtained by several types of annihilation 
final states which we will comment on later. 

If one uses the annihilation rate determined from relic DM calculation to calcu- 
late the excesses, one finds that the resulting excesses from the about propagation 
model calculations would be smaller by a factor of 100 to 1000. A large boost factor 
B is needed to explain the data. This is an important issue needs to be addressed 
before a consistent model can be constructed and tested further. 

This problem may be due our lacking of knowledge of DM substructures in 
our galaxy. Detailed calculations based on the N-body simulation shows that the 
boost factor from DM substructures is generally less than [134-137] ^ 10. Sev- 
eral other mechanisms have been proposed, including DM non-thermal production 
mechanism [131-133] the Sommerfeld effect [15] [25] [138-143] and the Breit-Wigner 
resonance enhancement effect [30] [144,145]. The non-thermal production mecha- 
nism is to detach the relic DM density from the annihilation rate producing the 
excesses. The annihilation rate is taken to be a parameter to be determined by the 

excess data. B = 1 can produce large enough excess. On the other hand, the 
Sommerfeld and Breit-Wigner mechanisms start with the annihilation rate deter- 
mined by relic DM density and dynamically determine the boost factor. We consider 
these latter two effects to be more natural. In the follow section we discuss these 
two mechanisms. 

3. The Sommerfeld and Breit-Wigner enhancement factors 
3.1. The Sommerfeld enhancement factor 

The Sommerfeld enhancement is a non-relativistic quantum mechanical effect [170]. 
Since at the epoch of relic DM got out of thermal equilibrium, the DM are non- 
relativistic, one can teat the problem with non-relativistic Schrodinger quantum 
mechanics. A large boost factor can be produced if DM interacts with a light par- 
ticle. 

The annihilation of DM is usually a short distance effect. Assuming that it hap- 
pens at the origin of a coordinate. The annihilate cross section can be approximated 
by a potential of the form Va5{r). Imagining now that the DM is moving in the z- 
direction, the wave-function can be written up to some overall normalization factor 
as, fplir) = e*'^^. It is obvious that the annihilation cross section is proportional to 
IV'fc(0)|^ since this factor represents the DM density at the origin where the anni- 
hilation occurs. The cross section cto due to the short distance interaction can be 
easily obtained by the usual scattering theoretical calculations. 

If in addition to the short distance interaction, there is an exchange of a massless 
particle between the two annihilation DMs, a long range interaction potential of the 
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form Vc = —a/r between the DM will be generated. Before the annihilation hap- 
pens, the 1/r long range force, as is well known, will distort the DM wave- function 
ip{r) and therefore modify the DM density at the origin. The cross section is then 
proportional to |^/'fe(0)p resulting in a modification factor, S = |^/'A;(0)P/|V'fc(0)P 
and the cross section is given by cr = (t"5. This factor is the Sommerfeld factor. 

The S factor due to a long range Coulomb like potential Vc{r) = —a/r can be 
obtained by solving the Schrodinger equation 

EMr) = (-^ V' +Vc{r))M^ , (H) 

where E — fc^/2/i > 0. /i = mDrriD / (m£) + mu) = tod/2 is the reduced mass in 
the center of mass frame of the two annihilating DM. 

This equation is a standard central force problem and can be solved by separation 
of variable in the form ?/>fe(r) = ^i"»i?;(r)r/"(6', </)) with Yl^{e,<p) being the 

spherical harmonic function. Ri{r) satisfies 

-2T'"t + - + ^^)) Ri{r) = . (12) 

dr dr / 

To obtain an analytic solution suitable for scattering problem, let us write 
Ri{r) = r'-e^^^ fi{r). The function fi{r) satisfies the confluent hypergeometric equa- 
tion [171] 

^^//(^) + m + 1) - z)-^Ji{z) -{1 + 1 + in)fi{z) = , (13) 

where z = —ikr and n = —fia/k. The solution is the regular confluent function 
fi{z) = iFi{l -|- 1 -|- in, 2{l + 1), z) since at r = 0, the wave function must be finite. 
Then 

Ri{r) = Cye'''\Fi{l + 1 + in,2{l + 1), -2ikr) , (14) 

To determine the constant C'l, one matches the asymptotic behavior at r — > oo 
of the partial wave decomposition of a plane wave scattered by a potential 



Mr,0) - e^"' + l—J2{2l + l)^-^Pi{cos0) 



oo 

— V(2/ + 1) (e"'+2'«' - e-^C^—'-] p,(cos 0) , (15) 



2^L ^ 



with the behavior of '(/'^(r, 9) — J^i CiRi{r) Pi {cos 6) at r — > oo 

gn7r/2p(2Z + 1) 1 



-i7r/2-nln(2fer)+r)j _ g-i(fcr-i7r/2-n ln(2fcr)+r)i 



^ ' {2kyT{l -f 1 -I- m) 2ikr V 

_ e"^/^r(2/ + 1) ,(i^/2+n In 2-m ^1 

'{2kyT{l + l + in) ' ^ ' 

^ / i(fcr-nln(A;r)-2nln2+2r;/) _ ^-i{kr-n\n{kr)~l7r) \ (^^qS ^) ' 
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where rji = argVil + 1 + in). 

Naively it seems that, keeping Ci independent of r, it is not possible to written 
the asymptotic form of eq.(16) into the corresponding coefficient (the same I) in eq. 
(fT5|) because the In(fcr) factor in the exponential. This is a well known fact which 
happens to the Coulomb potential scattering. The In(fcr) needs to be kept in the 
exponential due to the long range interaction nature of the potential. In this case 
the scattering phase 5i is equal to rji — nln2, and C; is then given by 

Ci = -^{2ik)^T{l + 1 + i„,)e""^/2e»("in2-m) (17) 
(zi). 

Note that whatever modification to the wave-function at r = 0, it only happens 
to the S-wave, I = case, since for I ^ 0, i?,(r) - r'e'''''iFi(Z + l + m, 2(? + l), -likr) 
is zero. One immediately finds that 

lAfelO) = r(l -f m)e-"'^/2gi(nln2-„o) (^8) 

Using the identity r(l -I- a;)r(l — x) — xtt/ sin(a;7r). one finally obtains 

a _ l^fc(0)P _ 2n7r ^ -qtt/i; 

|V'2(0)|2 e^"'^ - 1 e-"'^/" - 1 ■ ^ ' 

With the Sommerfeld factor included, the cross section is given by 

a{v)v = a''{v)v . (20) 

If a is positive, corresponding to an attractive potential, S is an enhancement factor. 
While for a negative a, corresponding to a repulsive potential, 5* is a suppression 
factor. Also it is velocity dependent. S plays the role of the boost factor. 

For S-wave annihilation, if the DM is not close to a resonant region, the cross 
section is almost a constant in DM velocity. Therefore boost factor for DM annihi- 
lation resulting from the relic DM is approximately the ratio of the DM velocity 
at the relic DM decoupling time in the early universe and the velocity Vp of DM in 
our galaxy halo at the present, that is i? ^ Vr/vp. The average Vr at the decoupling 
temperature Td is about ^j2Td/mD with Td/mo ~ 1/20 leading to Vr ~ 0.3. Model 
estimates show that Vp is about a few times 10"'*. With these numbers, one can see 
that the Sommerfeld enhancement can easily produce the needed large boost factor. 

In many practical situations, the interaction between DM is not mediated by 
a massless particle but a massive one, and the potential produced is a Yukawa 
potential Vyir) ~ — (Q;/r)e~™"''''. Here is the mass of the mediating particle. 
One needs to solve the differential equation in cq. p5|) with Vc [r) replaced by Vy (r) . 
Unfortunately with Vy (r) , it is not possible to obtain a simple analytic solution. To 
obtain the corresponding Sommerfeld enhancement factor, numerical calculation is 
needed. Fig. [3] shows the enhancement factor as a function of model parameters 
obtained in Ref. [25]. There are regions of parameters, resonant regions, where the 
enhancement factor can easily be as large as 1000. 

One can understand some qualitative features without detailed numerical calcu- 
lations. The Yukawa potential is not a long range interaction but one with a force 
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Fig. 3. The Sommerfeld factor contour as a function of = v/a and = m^/ amj:) . 

range cut-off at r ~ l/?Ti0 and therefore the enhancement factor will be smaller 
than that due to a Coulomb potential. For r < one can expand the expo- 

nential part of the potential in power of r. Keeping the leading correction, the 
potential is: VV(r) ~ —{a/r — am^). Comparing with eq. p^ . one notices that if 
ir^ » v1 = A\a\m^/mo, the correction term can be neglected and get back to 
the Coulomb potential case. At the velocity Vg the Sommerfeld enhancement starts 
to damp out. When the velocity slowed further down to where the corresponding 
deBrolie wave length l/moVd of the DM to be comparable with the force range 
the Sommerfeld enhancement saturates itself. This sets a natural limit of the 
enhancement [25], namu/m^. 

For a detailed discussion, an more accurate form for the averaged annihilation 
rate should be used which can be written as [172] 

(-') - -^1^ r (^) ds, (21) 

with 

a{s) = 2g^m^ ^ s - Am^ ■ av, (22) 

where x = nio /T and gi is the internal degrees of freedom of DM particle which is 
equal to 4 for a Dirac fermion. Ki{y) and i^2(y) are the modified Bessel functions 
of the second type. 

In the non-relativistic case, the above reduces to the Maxwell-Boltzmann aver- 
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age, 

(av) = J a{v)v^^v\-^'^''' , (23) 

where S is the average DM velocity. At the relic DM decoupling, it is y/2Td/mD, 
and at present in our galaxy halo, it is about 10^^. 

Finally to determine the parameters one should solve the standard Boltzmann 
equation governing the DM abundance [172] 



dY Xs{x) 2 -1^2 



<av> (y2 -Yio), (24) 



where Y = n/s{x) with n the DM number density and s{x) = 27r^(7*m|j/45.T'' the 
entropy density. H = ^ Ai:"^ j 45m|, jMp l is the Hubble constant evaluated at a; = 1. 
Here MpL ~ 1.22 x 10^^ GeV is the Planck mass, g^, is the total relativistic degrees 
of freedom. In the SM g^, = 106.75. The Y value at thermal equilibrium Yeq is 
given by: (45/4^2^7/2) (5,/g,)x3/2e-^ 

The rest challenge is more a particle physics one, finding a particle physics model 
which can have a particle (j) with small mass and interact with DM to produce the 
required boost factor. 

3.2. The Briet-Wigner resonant enhancement factor 

The Breit-Wigner enhancement mechanism can produce a large boost fac- 
tor 15,143,144 j£ ^Yie annihilation of DM is through exchange of a particle in the 
S-channel with mass close to 2m ij. Let us consider an example that the DM anni- 
hilates into lepton pairs by exchange of a S-channcl Z', ipip ^ Z' II through the 
following Lagrangian 

L={ag'i^^^'i, + g'hn)Zy (25) 
The interaction rate, to the leading order in velocity w, is given by [30] 

^(«)^' = - 7 ff'"'^2 . > (26) 

^(s-m|,)2 + r|,m|, 

where T z' is the total Z' decay width and s is the center of mass frame energy 
squared. 

If 2m £) is far away from mz' the interaction rate is almost a constant. After the 
parameters fixed by relic DM density requirement, the interaction rate is constrained 
to be too small to explain the data. A large boost factor is needed. The boost factor 
can arise from the fact that when the Z' mass is close to 2m u , the annihilation rate 
is close to the resonant point. In this case the interaction rate is very sensitive to 
velocity of the DM. To see this let us rewrite the above annihilation rate as 

16^ (^TW4FT^' ^ ^ 
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la' ^a' io* io= lo' 

x=my/T 



Fig. 4. The Breit-Wigner enhanced relative interaction rate {crv) / {av) ^=20 as a function of time 
X. here is mo. 

where we have used the non-relativistic Hmit of s = 4m|j + m|)W^, with 5 and 7 
defined as to|, = 4m|)(l — S), and 7^ = r|,(l — (5)/4m|). 

It is clear that for small enough S and 7, the annihilation rate is very sensi- 
tive to the velocity v. At lower velocity, the annihilation rate is enhanced. This 
results in a very different picture of DM annihilation than the case for the usual 
non-resonant annihilation where the annihilation rate is not sensitive to DM veloc- 
ity. The annihilation process does not freeze out even after the usual "freeze out" 
time in the non-resonant annihilation case due to the enhanced annihilation rate at 
lower energies in the early universe. To produce the observed DM relic density, the 
annihilation rate at zero temperature is required to be larger than the usual one, 
and therefor a boost factor. With appropriate 6 and 7, a large enough boost factor 
B can be produced. 

Once the interaction rate is obtained, one can use the formulae in eas. (|2ip and 
to obtain the boost factor B. A numerical evaluation obtained in Ref. [30] for 
the Breit-Wigner enhancement is shown in Fig. [4] 

4. Theoretical models for dark matter 
4.1. Dark matter models 

Many candidate WIMP DM models have been studied in the context of the recent 
e± data [25] [26-30] [31-48] [49-51] [52-54] [55-78] [79-95]. But not all of them 
can explain all the features of the excess in cosmic ray. It turns out that the 
spectrum shape of the observed excesses is easy to reproduce. For example, 
two boy annihilation directly into pair produces a hard spectrum which can be 
made consistent with ATIC data, but disfavored by FERMI-LAT and HESS data. 
If the two body final states are or pair, their subsequent decays into 
can produce a soft enough spectrum consistent with FERMI-LAT and HESS data. 
The large boost factor and excesses in energy spectrum but not in anti-proton 
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spectrum required from the data eliminate many candidate DM annihilation models. 
There are some doubts on the reliability of the null result for the anti-proton. But 
before further experimental results disprove this result, one cannot simple ignore it. 
In the following discussions, we will take these two as requirements for a successful 
model. 

The simplest DM model is the Darkon model [173]. This model is the simplest 
in the sense that it needs the least extension from the SM. It contains in addition 
to the SM particles, just one real SM singlet scalar. This singlet scalar plays the 
role of DM. The annihilation is through exchange of the SM Higgs boson. The mass 
can be as large as a TeV, but should not be too small (less than a GeV) in order 
not to produce too large detection cross section to be ruled out by data [174, 175] 
Although large boost factor can be induced if the Higgs mass is close to 2 times of 
the darkon mass through the Breit- Winger enhancement mechanism, it produces 
too many anti-protons in the cosmic ray after fitting the excess data. This is 
because that in this model the couplings of Higgs to SM fermions is proportional to 
fermion masses, the annihilation tends to favor heavy final states. Further extension 
of the model is needed to explain data. 

The most popular DM candidate is the lightest supersymmetric particle (LSP) 
in the minimal supersymmetric standard model (MSSM). The MSSM LSP has been 
proposed to explain the excesses [31-48] However the MSSM also has problem. 
The LSP is a linear combination of photino, zino and higgsino. It usually has a 
large hadronic annihilation fraction in conflict with the PAMELA data on anti- 
proton cosmic ray. Also in this model there is the problem to realize a large boost 
factor. Gravitino DM candidate also has similar problems. Extensions are needed. 
There are a few papers on this subject. The next MSSM (NMSSM) model has aU 
the right ingredients to explain the data. We will describe it more later. 

In universal extra dimension (UED) models, all SM particles live in extra dimen- 
sions. After compactification of the extra dimensions, there are Kaluza-Klein (KK) 
excitations generated. If the size the extra dimension R is of order 1/TeV. The KK 
mode can paly the role of DM to explain the excesses if a discrete symmetry 
is applied to make KK mode stable [52-54]. The KK DM can annihilate into SM 
particles. For example, the U{1)y gauge boson KK mode can annihilate through 
t-channel fermion KK mode into SM fermions, and produce excesses in cosmic 
ray [52-54]. However, they usually have a large hadronic fraction making the model 
troublesome with PAMELA anti-proton data. This problem can be remedied in the 
split-UED model where it is possible to split the lepton and quark KK masses such 
that the anti-proton cosmic ray is suppressed by larger quark KK mode masses [54] . 
However, these models all have problem to have a large boost factor. 

Many other models [25] [26-30] [55-78] [79-95] have been proposed to explain 
the data. There are basically two classes of models: a) kinematically limited light 
particle decay models [25] [55-78]; and b) leptophilic DM models [26-30] The light 
particle decay model a) requires the existence of a light particle with a mass less 
than the sum of proton and anti-proton masses and thus the light particle decays 
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predominantly into final states containing an e and/or a ^. If final states with 
such a particle is the dominant annihilation channel, the produced will be able 
to produce the excesses with appropriate mechanism to acquire a large boost 
factor. Exchange of a light enough particle between DM can produce a large boost 
factor through Sommerfeld enhancement mechanism. The Beit-Wigner mechanism 
can also supply the boost factor. The option b), leptophilic model, can be realized 
by interaction of DM being leptophilic (or hadrophobic) such that interactions of 
the particle mediating DM annihilate only have non-zero couplings to Icptons at tree 
level. In this case the mediating particle does not have to be light. Depending on 
how the annihilation occurs, it is also possible to have Sommerfeld or Brei-Wigner 
enhancement mechanism. In the following we discuss two models to illustrate how 
consistent models can be constructed. 



4.2. A leptophilic Z' model 

One of the following global symmetries in the SM can be gauged without gauge 
anomalies [176-180] 

a) hf. — Lfj, b) Lp — Lt, c) — . 

At the tree-level the Z' only couples to one of the pairs e and /i, e and r, and /i and 
T. If the Z' in one of these models is the mediating DM annihilation, the resulting 
final states are mainly leptonic states which can lead to excesses in observed in 
cosmic rays. Of course it requires that the Z' to couple to DM [14] [27] [30]. We 
assume that the DM field is a fermionic field ip with a non-trivial Li — Lj number 
a. To have an anomaly free theory, this DM field should have vector-like coupling 
to Z'. The Z' boson can develop a finite mass mz' from spontaneous U{l)Li-Lj 
symmetry breaking by a non-zero vacuum expectation values Vg of a scalar S with 
a non-trivial charge Li — Lj = b with m^, = b^g'^Vg. The Z' has the desired 
leptophilic couplings to fermions given by [30] 

L = -g'iay^-f^"^ + hj^h - Ijl^'lj + ua^Lv, - DjYLvj)Z'f, ■ (28) 

The relic DM density is controlled by annihilation of -ipij} — > Z'* — > IJi + ViVi. 
The interaction rate crw, with lepton masses neglected and summed over the two 
types of charged leptons and neutrinos, is given by 

_ 3 a'g'^ml 

^(.s-m|,)2 + r|,m|, ■ ^ > 

If the Z' mass is below the "0-0 threshold which we will assume, the dominant decay 
modes of Z' are Z' lili + DiVi, and Vz' is given by, neglecting lepton masses: 
Yz' = 3g'^mz'/l2Tr. 

In the above expressions for av and Tz' , it has been assumed that there are only 
left-handed light neutrinos. If there are light right-handed neutrinos, the factor 3 in 
these equations should be changed to 4. 
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When calculating relic DM density, one should use the above interaction rate. 
However, when calculating the spectrum, the above interaction rate needs to be 
multiplied by a factor of 2 /3 since the neutrinos in the final state do not contribute 
to the spectrum. 

The Breit-Wigner enhancement factor can be used to generate the needed large 
boost factor by requiring mz' to be close to 2m jj. The results obtained in Ref. [30] 
are shown in Fig. [5l The background is calculated using GALPROP package [3] 
with the diffusion 4- convection model parameters developed in Ref. [14]. Models 
a) and b), having hard in the final state, can fit ATIC and PAMELA data [30] 
with DM mass 1 TeV. Model c) can fit FERMLLAT data [30] with DM mass 1.5 
to 2 TeV. Once observational data finally settled down, some of the options can be 
further eliminated. 




Fig. 5. Left: positron fraction e+/(e+ -|- e~) predicted in the U(l)Li—Lj model compared with 
the observational data. Right: the normalized total electron spectrum of the model, compared with 
observational data, m/j = 1.5 TeV is used inthc figures. 



One may also construct similar model using Sommerfeld enhancement mecha- 
nism. A possible way of achieving this is to introduce a light real scaler cj) which car- 
ries no SM quantum numbers and therefore only couples to DM in the form a-ipipcj). 
Exchange this scalar will produce a Yukawa potential, — —{a/r)e~''"'>'^, be- 

tween DM. If the mass is small, a few GeV, there is enough parameter space 
where the model can fit the data. The results are similar to the case with Breit- 
Wigner enhancement mechanism. The field can in principle mix with the SM 
Higgs doublet through and then decay to SM particles. If the mixing is small 
enough it can avoid many potential phcnomenological constraints. 

The main difference for models with Breit-Wigner and Sommerfeld enhancement 
mechanisms is that in the former model niz' is fixed to be close to 2m]j, while the 
latter model niz' can have a much smaller mass leaving better chance for LHC study. 
Some of the discussions on LHC physics discussed in Ref. [27] can be applied. If there 
is kinetic mixing term [181] F^'^B^j^ between the field strength B^i, of the U{1)y 
and the field strength F'^'^ of the new U{l)Li-Lj , then there is chance to have more 
couplings to hadron states making LHC study more relevant. This scenario is worth 
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further investigation. 

Another interesting possibility is that the U{l)Li~Lj gauge boson Z' mass is 
actually very small [176-180] ( O(GeV)). With the Sommerfeld mechanism provides 
the boost factor, the is produced by t-channel — > Z' Z' first and then Z' 
decays into final product with . In this case the light Z' may be produced at low 
energy through [182] e+e^ Z'7. 

4.3. A light particle decay model 

As mentioned earlier that MSSM has problem accommodating excesses in cosmic 
ray data and there is the need to extend the model to explain data. It has been shown 
that the NMSSM [183-186] can have aU the required ingredients. The NMSSM has, 
in addition to the usual MSSM particle contents, a SM singlet chiral super-field 5. 
The fermionic partner of S is the singlino, x- The spin-0 partner S contains a scalar 
h = Re{S)/^/2 and a pseudoscalar a = /m(S')/V2. The allowed renormahzable 
super-potential Ws and soft SUSY breaking potential Vs are given by [49, 50] 

Ws = vis + ]^^isS^' + XHuHaS + ^kS^ , 

Vs = -[im^S-tS* + BsS^ + XAxHuHdS + kA^S^] + H.C. , (30) 

where Hu,d are the MSSM SM doublet super-fields. If a Z3 discrete symmetry is 
imposed on Wg and Vg, only 771^, A, k, A\ and A,^ are allowed. 

With suitable parameters in the model the singlino x can have a very small 
mixing with other SUSY neutral fermionic fields and is the lightest super particle 
(LSP) playing the role of DM. The singlino mainly annihilates into the scalar h 
and the pseudoscalar a, xx ~* through t and u channel exchanges of Xj and 
s-channel exchange of a. With non-zero XX^j XXo- and haa couplings, Th^x^ Taxx 
and guaa, respectively, the interaction rate is given by [49] 

a{xX ^ ah)v « [j^^alaaTaxx + Tlxx^axx ~ haaUxxT^xx^ (^l) 

The singlino x can also annihilate into hh and aa final state from t and u chan- 
nel X exchanges. But these contributions are suppressed by v'^ and do not play a 
substantial role. 

The h and a can mix with the Higgs and the pseudoscalar in the MSSM and 
therefore can decay into SM particles. If a has a mass iria below the threshold of 
hadronization, it will not decay into proton and anti-proton pair or even pions. If h 
mass JTift is larger than 2ma, h ^ aa is the main decay channel for h. This way the 
final products of xX a-h ^ aaa have no proton and anti-proton complying with 
PAMELA data. If rua is smaller than 2to^ but larger than 2m^, a predominately 
decays into fJ-^ i-i^ pairs and then fj. Vfj^ev^. In fact there is a hint from hyperCP 
measurement [187] that there is a light particle of mass 214 MeV from E+ pn~^ij,~ 
data and the NMSSM particle a fits that well [188]. The excess produced this way 
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will have soft spectra. It is not possible to produce a peak like spectrum in the 
ATIC data. But with x mass of order 2 TeV, the PAMEAL and FERMI-LAT 
spectra can be explained. 

Exchange of the scalar field h between singlino x produces an attractive Yukawa 
potential VY(r) — — (Q;/r)e~™''''. If rrih is light enough, less than a few GeV, a 
large Sommerfeld enhancement factor can be produced which then supply the much 
needed boost factor. This is a quit economic model. Fig.[6]shows how the PAMELA 
data can be fitted [49] with singlino mass of 600 GeV. But with x mass of order 
2TeV, although there may be the need of tuning the parameters, the PAMELA and 
FERMI-LAT spectra can be explained. 

If the mixing of a with the MSSM heavy pseudoscalar A is significant. The 
annihilation rate XX ^ A ^ ah can be large when the mass rriA of the field A is 
close to 2to^ . In this case the Breit-Wigner mechanism for the large boost factor B 
can be in operation, and also explain the data [51]. In this case is predicted to 
be 2mu which would in the range for 3 to 4 TeV. 

In this types of model, if the light intermediate state does couple to quarks, even 
it is kinematically limit to decay directly in to final states containing anti-proton 
or proton, there are off-shell contributions, but the rate is suppressed. 




Fig. 6. The cosmic ray positron fraction resulting from singlino DM annihilations. The e+ spec- 
trum is produced by xX ~* '"^ followed by /i —> aa and a fi'^ fi~ . mj, is taken to be 10 GeV. 
The dot-dashed line denotes the prediction from astrophysical secondary production alone. 



5. Discussions and Conclusions 

In previous discussions we have concentrated on explanation of excesses. Models 
constructed for this purpose have many other testable consequences. We list, with- 
out detailed discussions, some of the interesting subjects which can further reveal 
DM properties: a) Many of the models predict anti-proton excesses at higher en- 
ergies. Measurements of anti-proton with energy beyond the PAMELA range can 
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reveal more details of DM b) Almost all DM annihilation models predict certain 
level excesses of 7 ray. Further improved data on 7 ray can distinguish different 
models, c) Some of the models proposed, such as leptophlic models, predict very 
small cross sections for detection, but some other models with larger values. Direct 
detection of DM is certainly important to distinguish models, d) In some of the 
models there are light new particles. Negative results of direct detection of these 
particles can rule out some of these models, c) High energy collider search for DM 
or DM annihilation mediating particles or new light particles. LHC, ILC and CLIC 
can also play important role in distinguish different models. 

We conclude that DM annihilation can provide a consistent explanation for the 
recently observed excesses in cosmic ray. In order to cover the whole energy 
range of excesses observed, the DM mass must be as large as the highest energy 
observed showing excesses. The FERMI-LAT and HESS data then require the DM 
to be around 1.5 to 2 TeV. To produce large enough excesses with the annihilation 
mechanism, it requires modifications of the usual DM properties because that the 
same DM annihilation process is also required to produce the relic DM density 
in the early universe. With this constraint, a large boost factor in the range 100 
to 1000 is needed to explain data. This boost factor can be provided by particle 
physics effects. We have discussed two popular ones, the Sommerfeld and Breit- 
Wigner mechanisms. These two mechanisms have different consequences which can 
be distinguished by future experimental observations. The Brei-Wigner mechanism 
requires that the annihilation is through s-channel and the mediating particle has a 
mass close to two times of the DM mass. While the Sommerfeld mechanism requires 
the existence of a light particle of mass less than a few GeV. The PAMELA result 
of no anti-proton excesses further constrain how DM is annihilated. There are two 
classes of models, the leptophilic and the kinematically limited light particle decay 
models. The former requires the couplings cause the annihilation only have non-zero 
values for leptons, and the latter requires a new particle which the DM primarily 
annihilate into and this particle subsequently decays into leptons. Because this 
particle has a small mass which is kinematically forbidden to decay into proton 
or anti-proton and therefore explain the PAMELA null excess of anti-proton. All 
models modify what was called the usual WIMP DM. There are very different 
features for different types of models which can be tested. To further understand 
the properties of DM, more experimental observations are needed. 
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